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Abstract 

We present a new method for finding isolated exact solutions of a class of non-adiabatic Hamil- 
tonians of relevance to quantum òptics and allied areas. Central to our approach is the use of 
Bogoliubov transformations of the bosonic fields in the models. We demonstrate the simplicity 
and efnciency of this method by applying it to the Rabi Hamiltonian. 
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I Introduction 



There exists a class of simple, non-integrable, adiabatic Hamiltonians of the type that find application 
as models of light-matter interactions, for which it is possible to find exact isolated solutions. Generally 
these models involve some atòmic system, typically characterised by a simple two-level (or multi-level) 
system, interacting with a number of bosonic fields. Making the familiar rotating-wave approximation 
usually renders these models completcly soluble, but avoiding this approximation maintains the non- 
integrability of the models, and gives rise to the possibility of isolated exact solutions. This was 
first demonstrated for the Jahn-Teller model by Judd [Q , and these solutions are often referred to as 
Juddian solutions. Probably the simplest model for which these solutions have been found is the Rabi 
Hamiltonian (RH) , which describes a two-level atom interacting with a single-mode bosonic field via a 
dipole interaction Q . The Juddian solutions of the RH were first discovered by Reik and co-workers ||) , 
where they were seen to occur at the level crossings in the energy schema of the system. This turns 
out to be a general and important feature of these solutions. 

Apart from being of interest for what they teli us about the structure and symmetries of these models, 
the Juddian solutions are of considerable further value. Simple quantum òptics and related models, 
such as the RH, have long been utilised as test cases for various calculational techniques |Q-q|, and 
the possession of exact solutions facilitates their accurate assessment. Furthermore, the existence of 
isolated exact solutions in non-integrable quantum models is also of interest from the perspective of 
studying possible quantum chaos in such systems |?], || . In addition, it is hoped that these exact 
solutions may serve as useful starting points for perturbative treatments of the entire spectra of these 
models. 

In this paper we present a new and more general method for finding these isolated exact solutions, 
which we believe to have several advantages over the methods hitherto employed. Judd and Reik, 
working in the Bargmann representation, have used power series and Neumann series Ansàtze for the 
field mode. Neither of these approaches is particularly intuitive and the resulting àlgebra can become 
complicated. Küs and Lewenstein J^] have given a more concise approach which, as we describe later, 
is clearly related to the method that we describe here. For models such as the RH they used Bargmann 
representation Ansàtze for the field consisting of a finite number of bosonic excitations on top of a 
coherent state. They have also extended their method to some further systems, such as a three-level 
system and an auto-ionising ion. 

We believe that the method we outline in this paper is both more intuitive and more efficient that those 
discussed above, and that it reflects the essential physics of the systems to a greater degree. At the 
heart of the method is a simple canonical transformation of the bosonic field operators of the models. 
This transformation suggests the existence of exact solutions in a most direct manner. Our method 
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also has the advantage that it is easy to generalise, and is readily able to be extended to "two-photon" 
type interactions, in which two photons are required to induce an atòmic transition |ïo| . 
The remainder of this paper is organised as follows. In Sec. || we outline our method for finding the 
Juddian solutions. We describe in some detail the theory of Bogoliubov transformations of a boson 
mode and pay particular attention to their relation to the coherent and squeezed states. We then use 
these transformations to investigate the displaced and squeezed harmònic oscillators, to develop insight 



into the reasoning behind this approach. In Sec. [II we apply this method to the Rabi Hamiltonian, as 
an example of the use of this method. We then finish with some conclusions and indications of further 
work. 



II Methodology 

The models that we consider here consist of an atòmic system interacting with one or more bosonic 
modes. Each of these modes is described by annihilation and creation operators, b and respectively, 
which obey the usual commutation relation, 

M] = i. (i) 

In general the atòmic system will be described in terms of a set of matrices. For example, the two-level 
system in the RH is described by the SU(2) Pauli matrices. 

Our method for finding exact isolated solutions for such systems involves two components. First, one 
must choose an appropriate representation for the atòmic matrices and then, crucially, one performs 
a Bogoliubov transformation of the operators of the field mode. The nature of this transformation 
depends upon the type of interaction being considered and, with the correct choice of parameters, it 
leaves the Schròdinger equation in a form that admits exact solution with very simple Ansàtze. 

II. 1 Bogoliubov transformations 

A Bogoliubov transformation is a transformation from one description of a field mode in terms of the 
bosonic operators, b and 6', to a description in terms of new bosonic operators, b and say. This 
transformation is canonical so that the new operators obey the same commutation relation as the old 
ones, namely 

6,6+1 =1. ( 2 ) 
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The most general linear Bogoliubov transformation may be viewed as a rotation plus translation of 
the original oscillator Hilbert space to the new oscillator space, 

b=e- i0 {l-\a\ 2 y 1/2 (b-atf-z) 
6* = (1 - \<j\ 2 Y 1/2 (6t - a*b - z*) , (3) 

where a and z are complex numbers describing the amplitudes of the rotation and translation respec- 
tively. /3 is a simple, and usually rather unimportant, phase factor. From the outset it is important to 
note the restriction |er| < 1 in order to preserve the unitarity of the transformation. In the following 
we consider two specialisations of this transformation, namely a pure translation and a pure rotation. 
These transformations may be very simply related to the familiar coherent and squeezed states of 
quantum òptics and it is from this standpoint that we introduce the transformations. 

II. 2 Coherent bosons 

The usual Glauber coherent states, \z), may be defined as eigenkets of the single-mode bosonic anni- 
hilation operator |ïï| |, 

b\z) = z\z), (4) 
where z is a complex number. Such states are readily constructed as the following equivalent forms, 



1 |*| a _*6+ 



|0> (5) 
= e( 2bt -"* b )|0), (6) 

where we have normalised the coherent state such that (z\z) — 1. The exponential operator in Eq. (^) 
is denoted as follows, 

D(z) = e( z6+ - z * 6 ), (7) 

and is called the displacement operator. It is a unitary operator and we may readily use it to perform 
a unitary transformation of the field operators, 

D(z)bDÏ (z) = b-z = a 
L>(z)ò t L> t (z) = ò t -z* = a t (8) 

The operators D (z) form a representation of the Weyl (or Heisenberg-Weyl) group when multiplied by 
a trivial phase factor exp (i<f>), with cj> real. The operators a and obey the same commutator relation 
as the original operators, and thus we see this transformation to be a Bogoliubov transformation of 
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the type described as a pure translation above. Equations (Q) and (g) clearly imply 

a\z) = 0, (9) 

from which we see that the operator a annihilates the coherent state \z). Thus \z) may be considered 
as the vacuum state of the et- type bosons, and we rewrite it accordingly as |0; z) = \z), 

a\0;z} = 0. (10) 

We shall call these et-type bosons "coherent bosons" and write their number states as \n; z), such that 
a^a\n; z) — n\n; z). 

II. 3 Displaced harmònic oscillator 

The simplest application of the coherent bosons is to the displaced harmònic oscillator, 

Hv= l -(x + V2\) 2 + \p\ (11) 

in which the centre of the oscillator is shifted by an amount — y2A. Introducing the harmònic oscillator 
operators via 

p=-^(6t_ 6 ), (12) 

the Hamiltonian reads 

H D = & t ò + A(ò t +ò) + ^ + A 2 . (13) 

By performing a Bogoliubov transformation of the original bosonic operators to a new set of coherent 
bosons, and a, such that 

a = b + X, a = tf + \, (14) 

we may rewrite the Hamiltonian of Eq. ( |Ï3| ) in the form 

# D = aW^ (15) 

The eigenstates of this Hamiltonian are thus clearly seen to be the number states of the ct-type bosons, 
with corresponding 



5 



II. 4 Squeezed bosons 

Following Bishop and Vourdas [Q we construct the most general squeezed state, \z; p, 9, /3), by acting 
upon the bosonic vacuum |0) first with the displacement operator D (z) of Eq. (Q) and then with the 
pure squeezing operator 5 (p, 9, /?), 

\z;p,6,0) = S(p,e,p)D(z)\O), (16) 

The squeezing operator is given by 

5 (p, 9, f3) = exp ^-ipe-^òt 2 + \pe w b 2 ^j exp (i(3b%) , (17) 

where p,9,f3 are real parameters. It is a unitary operator, S* 5 — 1, and provides a representation of 
the group SU(1,1). Using a relationship given by Perelomov |Ï3[| , we are able to write the squeezing 
operator in the equivalent form 

S(a,(3) = exp Qafet 2 ) (l - \ a \'f b W exp (-^*6 2 ) exp (ifjtfb) , (18) 

where (3 is the same real parameter as above, and a is a complex number with modulus \o\ < 1, given 
by cr = — e~ 10 tanh (^p) • Using this expression, we can use the squeezing operator to make unitary 
transformations of the bosonic annihilation and creation operators, 

S(<r,0)bS* (a,0) = e~ (l - | cr | 2 ) _1/2 (fe - a tf) = c, 

S{<r,p)tf&(a,p) = e l/3 (l - |a| 2 )~ 1/2 (6 f - <r*6) = c f . (19) 

The operators c and satisfy the commutation relation [c, cH = 1 and thus the transformation 
fe, fe^ — > c, is a Bogoliubov transformation of the rotation type. From Eq. (|Ï9|), it follows that for 
any function of / (fe, 6*) 

Sf (b, 6^ & = f (c, ct) «- Sf (6, fet) = / (c, ct) 5. (20) 

Equation ( po|) implies that 56 = c5 and hence |z;cr/3) = \z;p9/3) are eigenstates of the annihilation 
operator c, 

c\z;af3) — cS(a,/3)\z} 
= S(a,(3)b\z) 

= z\z;af3). (21) 

If we consider the squeezed vacuum 5|0) = |0; cr,(3) = |0;er), we see that it is independent of (3 and 
that 

c|0;cr)=0. (22) 
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The number states of the c-type bosons are denoted \n; a/3), such that c^c|n; a/3) — n\n; a/3). We call 
the c-type bosons "squeezed" bosons. 

II. 5 Squeezed harmònic oscillator 

In position representation the squeezed harmònic oscillator has the form 

H s = \{l + 2X)x 2 + ^(l-2X)p 2 1 (23) 
where the real parameter A determines the degree of squeezing, with the restriction that |A| < i. 



Translating this into the Standard bosonic representation defined by Eq. (12) we have 



2 

We introduce squeezed c-type bosons defined by 



H s = tfb + \ + A f òt 2 + ò 2 ) . (24) 



and leave a real but undetermined for the moment. Making these substitutions into Eq. (|2j), we have 

H S = J]-^) { [-° + A + X ° 2 \ ( c2 + ct2 ) + + 1 - 4Aí7 ) ( ctc + l) } • ( 26 ) 
We eliminate the first term in this Hamiltonian by choosing 

-a + X + Xa 2 = 0, (27) 

giving, as one of the two solutions, 



(28) 



With this choice, the Hamiltonian becomes 



H s = {Jc+l (29) 



2 

The eigenstates of this Hamiltonian are clearly the number states of the squeezed c-type bosons, with 
eigenenergies 

E n = jn+ijíl (30) 

We note that the other solution of Eq. ( ]27| ) with a = l ea ds to the unphysical oscillator with 

Hg = — (c^c + i) íi, and since this Hamiltonian does not have square-integrable solutions, we discard 
it. 
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III Application to the Rabi Hamiltonian 



The Rabi Hamiltonian (RH) describes a two-level atom interacting with a single mode of quantised 
electromagnètic radiation via a dipole interaction |2fl . It is usually written in the form 



#Rabi = -W (T X + Wb% + g (V + b) (cr + + (7_) , 



(31) 



where coo is the atòmic level splitting, u is the frequency of the boson mode and g is the coupling 
strength of the atom to the held. The two-level atom is described by the Pauli pseudo-spin operators, 
which satisfy the SU(2) commutation relations 



[fffc) O'l] = 2Í£fcZmCT„ 



(32) 



where k,l,m € {x,y,z} with k ^ l and £kim is the anti-symmetric Levi-Civita symbol. We have 
defined the raising and lowering operators as 

o + = &x + i<Ty, o_=o x - ICJy. (33) 

It is convenient to rescale the Hamiltonian as iÏRabi = ^-^Rabi, where 

H Rahi = üia z + tfb + \(tf + b)a x , (34) 

and Q = ^ and A = = s -. There is a conserved parity l·l associated with the Hamiltonian. 

n 



i 



r,T \tfb+- (a z + i; 



exp 

— G z COS (7TÒ^6) 



(35) 



such that [i?Rabi,n] = 0. The parity operator l·l has two eigenvalues, 7r = ±1. The RH is not known 
to be integrable, but isolated exact solutions do exist. Here we use the technique outlined above to 
find these Juddian solutions. 

In order to do this we first require an appropriate matrix representation for the Pauli matrices, which 
for this model is one in which a x is diagonal. We shall use 



(36) 



In terms of the two-component wavefunction, I*) — the time-independent Schròdinger equation 

for the system, iÍR a bi|í') = E\^>), then reads 





1 







i 







1 























-1 


-i 




1 






(1-1*2) + (ò t ò+A(ò t +ò) = 

cDl^!) + (6+6 — A (6 f + ò) - £7) |* 2 > = 0. 



(37) 
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We now make the Bogoliubov transformation to the coherent bosons, and et, specified by 

a t = ò f - A, a = b — X. (38) 

The vacuum state of these bosons is the coherent state |A). It should be noted that this choice of 
transformation may be intuited from considering the w = limit of the Hamiltonian, where the same 
transformation is used to solve the model exactly in this limit, which is essentially equivalent to the 
displaced oscillator considered earlier. With this transformation Eqs. ( |37| ) become 

w|* 2 ) + {a f a + 2A (a f + a) + 3A 2 - E) = 

Q^x) + {a)a- \ 2 - E}\^> 2 ) =0, (39) 

where the kets 1^1,2) are now in the transformed representation. For these kets we choose the Ansatz 

N-l N-l i j\n 

- y>>;AH Vp„-^i r |0;A>-P JV _ 1 (at)|0;A); 

n=0 n=0 v 

1*2) = Vç„|n;A) = Vç I1 ^i r |0;A)=Q J v(at)|0;A); (40) 

n=0 n=0 v 

where \n; A) are number states of the coherent bosons, a^a|n; A) = n\n; A), and we have introduced the 
polynomials Pn-i and Qm of order N — 1 and N respectively. Making these substitutions we have 

N N-l 



U) 

n=0 n=0 

N-l N-l 



Y,Qn\n;X) + ^p n {n + 2,X 2 -E) |n;A) 

n=0 
N-l 

2A PnVn+l\n + 1; A) + 2A ^ p n Vn\n - 1; A) = 0, 

n—O n—1 

N-l N 

w^p„|n;A) + In (n - X 2 - E) \n\ X) = 0. (41) 

n=0 n=0 



Considering the highest number state, \N\ A), in the second of these equations, we see that for this 
equation to hold we requirc 

(N - A 2 - E) q N = 0. (42) 

Since c/tv 7^ by Ansatz, we obtain a determination of the energy 

E = N-X 2 . (43) 

This equation identifies the Juddian baseline energies, along which the Juddian solutions lie. Com- 
paring the coefficients of the remaining number states gives us 2N + 1 linear equations for the 2N + 1 
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cocfficients (p m ,0 < m < N ~ 1) and (<?fc,0 < k < N). To obtain non-trivial solutions, we clearly re- 
quire the determinant of this equation set to be zero. This gives the compatibility condition, providing 
the locations of the Juddian points. The fïrst two conditions (N = 1, 2) have the explícit forms 

w 2 +4A 2 = l, foriV = l, (44) 

+ (12A 2 - 5) Ço 2 + 32A 4 - 32A 2 + 4 = 0, for N = 2, (45) 

as have been given by Küs and Lewenstein Q. Thus, for a given N, we have a polynomial of iVth 
order in A 2 and Cj 2 . Each of these has N roots for A 2 in terms of w 2 , which all turn out to be real, thus 
giving the location of N Juddian solutions. Bcfore we look at these results, it is of interest to consider 
the other possible type of finite Ansatz at the Juddian points. These are found by using the coherent 
bosons 

a f = + A, a = b + A, (46) 
and interchanging the roles of and l^)- 



III. 1 Results 

By solving the complementary conditions we have calculated the first ten Juddian points for the 
resonant RH. These are displayed in Table |ï|, listed to 10 decimal places. 

The location of these Juddian points in the energy schema of the Hamiltonian is displayed in Figure [ï], 
where the schema was obtained by approximate numerical diagonalisation via a Standard configuration- 
interaction method, using a basis size of the lowest 101 harmònic oscillator states [Q. Also plotted 
are the Juddian baselines from Eq. (|4^). 

From Fig. [ï] we see that the Juddian points occur at the level crossings in this diagram. Thus we 
see that they occur when two solutions of different parity II become degenerate in energy, and this 
degeneracy is the key to the existence of the Juddian solutions. The coherent-boson number states 
\n; A) are not eigenstates of II, and thus the Ansatz ( fïo|) is not of definite parity. It is precisely because 
we can construct wavefunctions of mixed parity that allows us to find such simple Ansàtze at the 
Juddian points. 

We are now able to make explicit the connection between this method and that used by Küs and 
Lewenstein Q in investigating the RH. They worked in the Bargmann representation ]ïq| , in which 
the bosonic operators are represented by 

òt^z; b^^. (47) 
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and postulated the following forms for the two components of the wavefunction: 

* 1 (z) = e-l z l 2 / 2 (z|vl/ 1 ) =P N _ 1 (z)e Az , 
* 2 (z) = e -W 2 / 2 (z|v]/ 2 ) = Q N (z) e Az , (48) 

where P/v-i (z) and Qn (z) are polynomials in z of order N — 1 and N respectively. Bearing in mind 
the form of the coherent state (^), these wavefunctions are simply seen to be of the form of polynomials 
in the bosonic creation operator, , acting upon a coherent state of amplitude A. In our Ansatz (|40|). 
we have the same coherent state but now being acted upon by polynomials in = (6^ — A) , which 
shares a closer connection to the coherent state than . 

The polynomials of Küs and Lewenstein, are simply related to those of Ansatz (^) by P/v (z) = 
P/v (z — A). In the present case where we have only used displacemcnts of the boson mode, the 
difference between the two approaches is thus minimal. However, this is not the case when we require 
the use of squeezed bosons. Generally, an Ansatz posited in the squeezed representation would contain 
polynomials of the form P/v (c'j, where c 1 is the creation operator of the squeezed bosons. The 
analogous Ansatz to Eq. (^) would still contain a polynomial in z, pv (z) say If we assume the 
simplest type of squeezing and write c< = (b^ + <rò) / (VT— c 2 ) as in Eq. (p5|), then the Küs and 
Lewenstein polynomial can be written 

P N (z) = Pv [y^j , (49) 

which, crucially, contains both z and its derivative, and although formal relationship do exist between 
the polynomials of the two methods, these relationships are generally not trivial, especially if one con- 
siders the more general form of the Bogoliubov transformation. So the Ansàtze of the two methods are 
seen to be significantly different, and we conjecture that the one described here has several advantages 
which we shall discuss in the conclusion. 



IV Conclusions 

We have presented a method for finding isolated exact solutions of a class of non-adiabatic models, of 
the type frequently used in quantum òptics and related fields. 

Compared with the original approaches of Judd and Reik, the above method is more transparent and 
considerably simpler, advantages that it shares with the technique of Küs and Lewenstein. However, 
we believe that the use of transformed bosons is more obviously physically meaningful than the use of 
wavefunctions in Bargmann space, especially given the connection of these bosons to the coherent and 
squeezed states, so important in quantum òptics. 
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As an example of the use of this technique, we have applied it to the Rabi Hamiltonian and obtained in 
a simple fashion the known Juddian solutions of this model. In this example, we have used the coherent 
bosons to obtain Juddian solutions for a problem with an interaction of the type A (b^ + ò) a x . It is 
hopefully now clear how one may apply this method to further problems containing the same type 
of interaction. We have not as yet mentioned the application of the squeezed bosons in performing 
this kind of calculation. This second type of Bogoliubov transformation is useful in finding Juddian 
solutions of models containing two-photon type interactions. An obvious example is the two-photon 
Rabi Hamiltonian |Ï6| ] , which has the Hamiltonian 

H = ü<r z + bïb + \(bï 2 + b 2 ^a x . (50) 

Using squeezed bosons we are able to obtain a set of Juddian solutions for this model and these results 
will be discussed in a future publication. 

Due to the intuitive nature and simplicity of this technique it is easy to extend to other systems. For 
example, in view of their mode of construction we expect that our displaced and squeezed coherent 
states will be of particular use in any quantum field theory that has underlying dynamical symmetry 
of the Weyl group or the SU(1,1) group, or to which the (inhomogeneous or homogeneous) Bogoliubov 
transformation may be profitably applied. The obvious group-theoretical foundations of the technique 
also point the way to other approximations, since, for example in the squeezed (two-photon) case, 
SU(1,1) is not the only relevant group. Thus, the three-dimensional Lorentz group SO(2,l), which is 
the group of rotations in three-dimensional Minkowski space with two space and one time dimensions, 
is locally isomorphic to SU(1,1). Similarly, both the groups SL(2,R) of real second order matrices with 
unit determinant and the symplectic group Sp(2,R) are also locally isomorphic to SU(1,1). 
One may also readily generalise the current approach for the two-level models involving linear or 
quadratic interactions with a single boson (or canonical quantum mode) to the corresponding case 
of linear or bilinear interactions involving several distinct bosons or modes. For the linear models 
involving only displacements this is essentially trivial. However, for models involving squeezing, in the 
case of n bosons or modes the various bilinear produets of operators b\ bj , bi bj and b\bj , i, j = 1, 2, . . . , n 
now form a realisation of the highcr symplectic àlgebra Sp(2n, M). As before one can simply construct 
a unitary representation of this group by exponentiating the skew-adjoint operators in the àlgebra. 
For example, Bishop and Vourdas |TÏ]] have shown explicitly how to construct the most general two- 
mode squeezed states associated with a unitary representation of the group Sp(4,R). Once again 
such states are the ordinary coherent states with respect to the new destruction operators c\ and C2, 
which are themselves general linear Bogoliubov transformations of the original destruction operators 
b%, 62 and their Hermitian-conjugate creation operators b\, b\. The Sp(4,K) àlgebra has various 
subalgebras corresponding to different sorts of linear pairing terms. For example, whereas the single- 
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mode paring operators = \ (b\j ; K^} = \b\\ = \b\bi + \ for i = 1,2 correspond to 

the so-called §) representations of SU(1,1), the mixed pairing operators L + = b\b 2 ; L_ = Ò1Ò2; 
Lq = i ^&í&i + &2^2 + 1^ correspond to the discrete-series representation of SU(1,1). By contrast, the 
mixed pairing operators J+ = b\b 2 ; J- = b\b\; Jq = | (b\bi — b\b2^j correspond to the (Schwinger 
representation of) the angular momentum subalgebra SU(2). 

Bishop and Vourdas have shown in a separate publication ]ï^] how squeezed (pair) coherent states 
can also be used in connection with a rather broad class of quantum Lagrangians which include the 
dampcd harmònic oscillator, and hence with problems involving "quantum friction" or fluctuation- 
dissipation phenomena in general. Within quantum òptics for example, the quantum theory of làsers 
and photon detection provide obvious applications. Such problems can now also usefully be extended 
by our present treatment to the case of such damped systems coupled to two level àtoms. 
The possibility of using these solutions as the basis of a perturbative approach extends the method away 
from just the isolated exact points to the remainder of the spectrum of the system. The properties of 
such an approach are yet to be investigated. Finally we note that the extension to similar single-mode 
or multi-mode systems as considered above coupled to n-level àtoms with n > 2 is also straightforward 
in principio. 
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9 


E 


N 


0.2165063510 


0.8125000000 


1 


0.1661640732 


1.8895580031 


2 


0.4460403578 


1.2041919969 


2 


0.1400889590 


2.9215003343 


3 


0.3664714887 


2.4627945920 


3 


0.6163829153 


1.4802884071 


3 


0.1234229399 


3.9390671161 


4 


0.3199075781 


3.5906365658 


4 


0.5243395120 


2.9002723045 


4 


0.7582492415 


1.7002323511 


4 



Table 1: The couplings, energies, and N, of the first ten Juddian points of the resonant Rabi Hamil- 
tonian (u = uç, = 1). 
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g 

Figure 1: The first ten Juddian points of the Rabi Hamiltonian (diamonds). Also plotted are the 
energy levels obtained by numerical diagonalisation (dark lines), and the Juddian base-lines (light 
lines). The Hamiltonian is resonant; to = Wq = 1. 
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